Relaxation dynamics of the electron distribution in the Coulomb blockade problem 
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We study the relaxation dynamics of electron distribution function on the island of a single electron 
transistor. We focus on the regime of not very low temperatures in which an electron coherence 
can be neglected but quantum fluctuations of charge are strong due to Coulomb interaction. The 
quantum kinetic equation governing evolution of the electron distribution function due to escape 
of electrons to the reservoirs is derived. Analytical solutions for time-dependence of the electron 
distribution are obtained in the regimes of weak and strong Coulomb blockade. We find that usual 
exponential in time relaxation is strongly modified due to the presence of Coulomb interaction. 
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I. INTRODUCTION 

Recently, the phenomenon of Coulomb blockade in sin- 
gle electron devices^— has come into the focus of the 
field of thermoelectricity.— Among significant experimen- 
tal achievements one can list development of the Coulomb 
blockade thermometer^ the thermal rectifier on the basis 
of a quantum dotj^ new technique to measure tempera- 
ture gradients across a quantum dot^ etc. The standard 
characteristic of thermoelectric performance is the fig- 
ure of merit which involves the product of conductance, 
thcrmopower squared and inverse thermal conductance. 
Measurements of thermopower and thermal conductance 
in single electron transistors (SET) and quantum dots 
have been performed during the last decade at different 
temperature regimes^— The theory of the thermoelec- 
tric effects in electron devices has been put forward in 
Refs. (j^fl3j . During the last decade the thermopower 
and thermal conductance have been studied in single 
electron transistors and quantum dotSfi 4 ^— and in gran- 
ular metals2i~— in various regimes. However, the ther- 
mopower and thermal conductance arc linear response 
parameters and, therefore, describe the equilibrium prop- 
erties of a system only. 

In contrast, our work is focused on properties of single 
electron devices in the out-of-equilibrium regime which 
has attracted a lot of theoretical interest recently. In 
particular, the conductance of a quantum dot under ac 
pumping in the stationary non-equilibrium state was ob- 
tained in Ref . [2fj| | , the current noise of an ac- biased quan- 
tum dot was studied in Ref. [27j . the non-equilibrium 
dephasing rate and zero-bias anomaly in single electron 
transistor was computed in Ref. [28[, the statistics of 
temperature and current fluctuations in the fully out- 
of equilibrium single electron transistor was investigated 
in Refs. 29ll30j . and the extension of the P (E)-theoiy^- 
to the out-of-equilibrium regime has been developed in 
Refs. [32II33I ]. However, these works dealt with regimes 
when the distribution function of electrons in a quantum 



dot or an island of single electron transistor were fixed 
by external sources, e.g., ac or dc bias voltage. 

In this paper we address a different question: how does 
an electron distribution function once being prepared re- 
laxes toward the equilibrium state in the Coulomb block- 
ade problem. Apart from general physical interest in un- 
derstanding of a non-equilibrium regime, the answer to 
this question is important for the field of electron ther- 
mometry^ 

We consider the simplest system: single electron tran- 
sistor. The set-up is shown schematically in Fig. [T] 
Metallic island is coupled to an equilibrium electron 
reservoirs via tunneling junctions. Depending on the 
task, the reservoirs and the island may be kept at differ- 
ent temperatures (T\, T r respectively), different chemical 
potentials (constant or varying in time: /i;(t), ^ r (t)) or 
quasi-stationary gate voltage (U g (t) = Uq + U u cosuit) 
may be applied to the system. The physics of the sys- 
tem is governed by several energy scales: the Thoulcss 
energy of an island E^h, the charging energy E c , and the 
mean single-particle level spacing 8. Throughout the pa- 
per the Thouless energy is considered to be the largest 
scale in the problem. This allows us to treat the metallic 
island as a zero dimensional object with vanishing in- 
ternal resistance. The dimcnsionless total conductance 
(in units e 2 /h) of tunneling junctions g is an essential 
control parameter. The junctions are assumed to have 
a large number of channels but the conductance of each 
one is assumed to be small gip <C 1. The temperature 
is assumed to be low enough: T -C max{l, g}E c in order 
to keep electrons strongly correlated due to Coulomb in- 
teraction. At low temperatures the interplay of Coulomb 
interaction and electron coherence dominates the physics 
of single electron devices. Account of both effects is a 
formidable undertaking. Therefore, we are going to re- 
strict ourselves to the regime of not very low tempera- 
tures in which an electron coherence can be neglected 
but quantum fluctuations of charge are strong due to 
Coulomb interaction. Then, the physics of the system is 
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FIG. 1: The sketch of a SET device. The leads are kept 
at different temperatures (chemical potentials) inducing heat 
(electric) currents. 



adequately described in the framework of Ambcgaokar- 
Eckern-Schon (AES) effective action^ The effective ac- 
tion describes the system in terms of bosonic field <p, 
which is usually termed as the plasmon field. Its time 
derivative is interpreted as fluctuating electric potential 
of electrons inside the island. AES-approach has well- 
known limitations. Deriving AES action one assumes 
that the products of electron Green's functions averaged 
over disorder are substituted with products of disorder- 
averaged Green's functions in every calculation. That is 
why the processes of phase-coherent multiple impurity 
scattering inside the island are left out. The limitations 
in the re gim e g > 1 and j< 1 were discussed in detail 
in Refs. 35 1 and (36|, respectively. It was shown that 
at temperatures T 3> maxjl, g}5, AES-action approach 
is justified. Following Ref. [33], we shall term the tem- 
perature range, maxjl, g}E c T maxjl, g}8 as an 
interaction without coherence regime. This 'interaction 
without coherence' regime is an attainable experimental 
reality, e.g. in experiments reported in Refs. Q and [To[ 
the necesary conditions were satisfied. 

In the case of strong Coulomb blockade (g <C 1), the 
theoretical study of relaxation of an electron distribu- 
tion in the interaction without coherence regime has been 
done before for a single quantum dot [38| and for an ID 
array of quantum dots [39(. However, the considerations 
of Ref. [38[ have been restricted by assumptions that i) 
the electron distribution is the Fermi function with some 
temperature different from the equilibrium one; ii) trans- 
port is dominated by co-tunneling processes (Coulomb 
valley regime); iii) temperatures of the island and the 
reservoirs are close to each other. 

In the present paper, we undertake the analysis of re- 
laxation of an electron distribution function which is free 
of above-mentioned restrictions. 

Since we are going to capture non-equlibrium physics, 
we employ the formalism of AES-action in its out-of- 
cquilibrium form throughout the paper. We supplement 
it with quantum kinetic equation to explore relaxation 
dynamics of electron distribution. For a SET with large 
number of tunneling channels we derive the quantum ki- 
netic equation with the collision integral due to escape of 
electrons to the reservoirs. It is valid in the entire s pan 
of values of g and generalizes the one obtained in Ref. [26[ 



for sequential tunneling (first order in g) and cotunneling 
(second order in g) approximations in the framework of 
the orthodox theory of the Coulomb blockade. In fact, 
our collision integral is always an infinite series in pow- 
ers of g. Indeed, each tunneling event is accompanied 
by the radiation of a plasmon. That is why the collision 
integral becomes of the infinite order in the distribution 
function of electrons inside the island. This situation is 
entirely different from the one in Fermi liquid and leads 
to non-trivial relaxation. 

As a test of the quantum kinetic equation, in the 
regime of linear response we derive analytical expressions 
for transport coefficients: conductance, thermal conduc- 
tance and the response of electric current to tempera- 
ture difference. In the regime of weak Coulomb blockade 
(g ^> 1) we establish the following new results for the 
transport coefficients: i) the conductance and thermal 
conductance violate Wiedemann-Franz law, and devia- 
tion of the Lorentz ratio C from value 7r 2 /3e 2 demon- 
strates weak periodic dependence on the gate voltage; 
ii) the thcrmopower weakly oscillates with the gate volt- 
age around zero value. Weak oscillations of the Lorentz 
ratio and thermopower with the gate voltage found in 
the regime g 1 are manifestation of the known gate- 
voltage dependence of these quantities^— in the strong 
Coulomb blockade regime, g <C 1. 

In weak and strong Coulomb blockade regimes we have 
employed the quantum kinetic equation to solve the re- 
laxation of the electron distribution in two cases: i) the 
distribution of electrons inside the island is the Fermi- 
function with some temperature; ii) the distribution func- 
tion of electrons inside the island is arbitrary. In the for- 
mer case we have managed to extract the relaxation dy- 
namics of the electron temperature; in the latter case we 
have obtained evolution of a distribution function itself. 
In both cases we assumed that electron escape to reser- 
voirs is the primary relaxation mechanism. In general, 
the collision integral in the quantum kinetic equation is 
non-local in energy due to inelastic nature of tunneling 
processes: the radiation of plasmon always accompanies 
the tunneling event. In a number of wide parametric 
regimes: weak Coulomb blockade and Coulomb peak in 
the strong Coulomb blockade, the kernel of the quantum 
kinetic equation acquires a quasi-elastic form. However, 
the collision integral remains non-local in energy due to 
rcnormalization effects in these cases. The co-tunneling 
regime is qualitatively different: the kernel of the colli- 
sion integral is entirely inelastic. 

Our new result is that despite quasi-elastic form of the 
collision integral, strong Coulomb interaction dramati- 
cally changes the relaxation laws comparing to simple 
exponential ones expected from golden-rule type argu- 
ments. They suggest that electron relaxation rate is to 
be proportional to the width of electrons' levels inside the 
island, g5, prompting simple exponential relaxation. The 
rcnormalization effects due to Coulomb interaction make 
the width of electrons' levels dependent on the electron 
distribution and lead to the non-exponential relaxation 



3 



laws. For example, in the regime of the sequential tun- 
neling, we have discovered that there is a time regime in 
which relaxation of the electron temperature in a SET 
island is independent of the tunneling conductance g. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the hamiltonian and essential parameters of the 
problem. Sec. IIIII is devoted to the out-of-equilibrium 
AES-model and to derivation of the quantum kinetic 
equation. Sec. [TV] is devoted to derivation of general 
expressions for the linear response coefficients. The re- 
laxation dynamics of electrons in the island is explored in 
the weak (g 3> 1) and strong (g <C 1) Coulomb blockade 
regimes in Sec. IVllVIl Discussion of the results, compar- 
ision with other relaxation mechanisms, different from 
electron escape to reservoirs and conclusions are pre- 
sented in Sec. IVIII 



II. FORMALISM 

A SET is described by the Hamiltonian 

H = H +H c + H t , (1) 

where 



(») Wt (i) 



(2) 



describes free electrons in the leads and the island, H c de- 
scribes Coulomb interaction of carriers in the island, and 
H t describes the tunneling. Here operators a£ (^1) cre ~ 
ate a carrier in the i-th lead (island). Then, the tunneling 
hamiltonian is 



ff t = E€4 i)t ^+H.c. 



(3) 



The charging Hamiltonian of electrons in the box is taken 
in the capacitive form: 



H c = E c (n d - q) 



(4) 



The energies e% , are accounted from the Fermi level, 
and the delta-functions should be smoothed on the scale 
5E, such that max{6, <5 (i ' r) } < 5E < T. Here, S and 
S^ l,r ^ stand for mean level spacing of single-particle states 
on the island and reservoirs, respectively. The classical 
dimensionless conductance (in units e 2 /h) of the junction 
between a reservoir and the island can be expressed as 
follows^ 



91 +9r, 



9l,r = Y. g kk ] =^2- 



act 



(8) 



Therefore, each non-zero eigenvalue of or g^ corre- 
sponds to the transmittance of some 'transport' channel 
between a reservoir and the island^ The effective num- 
ber of these 'transport' channels (N^) is given by 

2 



N. 



« 



l^9kk 

k 



/ , yaa 



ch 



9kk'9k'k 2^, 9aa'9 a 'a 
kk' olol' 

,(') 



(9) 



The effective dimensionless conductance g^J of a 'trans- 
port' channel can be written as follows 



l^9kk'9k'k 9 aa '9a'a 



(i) kk 

9ch = — 



T,9kk 

k 



Eg 



(10) 



The dimensionless conductance g then becomes 

9 = 9l+9r, gi.r - .y ch - ch 
In what follows we will always assume 



g^ r) N^ r) . (11) 



9% « 1- 



JV« » 1. 



(12) 



Notice that under these circumstances the conductances 
gi^ r can still be large provided the effective number of 

channels N^t r ' is sufficiently large. 



Here E c = e 2 /(2C) denotes the charging energy, and fid 
is an operator of a particle number in the island: 



(5) 



To characterize the tunneling it is convenient to intro- 
duce the following hermitean matrices: 

til = (2nf [6{e^)8(e^)] 1/2 ^t^S^)t^ (6) 

CX. 

til = (2n? teeW)^)]^!;*^?)*^, (7) 



III. ACTION AND KINETIC EQUAITONS 

A. AES-action 

To tackle the system which is out of equilibrium we 
have to employ essentially non-equilibrium formalism. 
Keldysh technique is thus the only way through. We em- 
ploy Keldysh form of AES-action (we sketch the known 
details of derivation in Appendix O 



S — S c -+- Sd 



(13) 



where 



The first of them acting in the Hilbert space of the states 
of the lead, the second - in the space of the islands states. 



S c = — J tpc'Pqdt + 2q J ip q dt. 



(14) 
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Here ip CtQ = ((p+ ± (fi-)/2 with <p± denoting bosonic field 
on both branches of Keldysh contour. Physically, the 
bosonic field is associated with the fluctuating electric 
potential on the island. In terms of classic and quantum 
boson exponents 



the dissipative part of AES-action reads: 



(15) 



S d 



(x c (t)x q (t)) 







n A (t,t') 
n K (t,t') 



Xc ^ \ dm' 
x q (t'); 



(16) 



Here H R > A ' K are corresponding components of electron 
polarization operator in the Keldysh space. They are 
given by a standard formulae presented for reference in 
Appendix [21 In a case of constant density of states 
(DOS) in the island and leads the kernel of the AES 
action can be simplified: 



n 



R.A.K 



(t,f) 



— n 

2tt 



R.A.K 



(r)e 



— iuj{t — t ! ) 



H^ A (r)=TiJ2 9 j J [F £ d (r) - F?_ u (t)} 

9a 



nff(r) 



^ 9 



2tt' 
ds 



{1-F?(t)F?_ u (t)) — . 

Zir 



(17) 
(18) 

(19) 



Here we define a slow time t = (t + t')/2. Function F £ (t) 
is given in terms of the Wigner transform f £ (r) of the 
electron distribution function f(t, t 1 ): F 6 (t) = 1 — 2f e (r). 
H R ' A ' K (t) arc Wigner transforms of corresponding func- 
tions in time domain. 

As seen from the structure of the r.h.s. of (TT91 . it is 
suitable to introduce function FJ = (J2 a 9aF")/g which 
may be called the effective distribution function of two 
reservoirs. It is that combination of reservoir distribution 
functions that enters all the following equations of the 
paper. 

Although Eq. (TT51) is exact we neglect all derivatives 
with respect to slow time r in Eq. (TIT)]) . It is also conve- 
nient to introduce function B u (r) in accordance with 



n^(r)=2iIm^(T)B u (r) 



(20) 



relating Keldysh and retarded (advanced) components of 
polarization operator. The function B^ir) plays a role of 
a distributiion function for electron-hole excitations. In 
the equilibrium it is given by coth(w/2T). 

In what follows we assume that electrons in the leads 
are locally thermalized such that f\' r (r) arc the Fcrmi- 
functions. Depending on the parameters of the model, 
both quasi-equilibrium and non-equilibrium regimes can 
exist in the island. Therefore, we assume F d {r) to be an 
arbitrary function slow varying with time r, and derive 
the kinetic equation for F d (r) by which the AES action 
should be supplemented. 






FIG. 2: Feynman diagrams for fermion self-energy: a) second 
order in Ht; b) fourth order in Ht- 



B. Kinetic equations 

The starting point for deriving kinetic equation for a 
SET is the Dyson equation for the Keldysh component 
of electron's Green's function^ 



(d t + d t ,)F d (t,t') 



Air 2 



Z R ■ F d + F d ■ 



t,t' 
(21) 

Here, va = S^ 1 = J2 a ^( £ a) ^ s an averaged single parti- 
cle density of states in the island and Yi K ' R ' A are the 
components of self-energy in Keldysh space. To the 
second-order in tunneling Hamiltonian Ht (lowest order 
in l/iV c h) the Wigner transform of the self-energy shown 
in Fig. [5^ reads (see Appendix [B]) 



K(r)=ng 



2 J 2tt 

duj 



R,Ai 



2tt 



(22) 



Here we perform Wigner transform of the exact sel- 
energies and introduce the correlation functions of boson 
exponents: 



V R (t,t') 
V A {t,t') 
V K {t,t') 



-i(X c (t)X q (t')), 
-i(X q {t)X c (t')}, 
-i{X c (t)X c {t')}. 



(23) 



It is convenient to parametrize T) K (t,t') via the boson 
distribution function B(t,t'): 



Df ( , = [V H -B-B-V J 



(24) 



It is worthwhile to mention that the next ( fourth order 
in Ht) contribution to the self-energy which is shown in 
Fig- lib is of the order g 2 /N c h- This correction to the self- 
energy is of the same order as terms omitted in the course 
of derivation of the AES action ([TB")) . In the considered 
limit N c h ^> 1 it can be safely neglected. 

Performing Wigner transform of Eq. (j2"Tj) and neglect- 
ing all slow time derivatives in its r.h.s. we obtain the 
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quantum kinetic equation for the distribution function of 
electrons on the island of a SET: 



a—l,r 

X {(Ff- U (r) - F*{t))b u {t) + 1 - F £ d (T)F*_Jr)}. 

(25) 

This quantum kinetic equation constitutes one of the 
main results of the present paper. It describes evolu- 
tion of the distribution function F£(t) of electrons in the 
island due to interaction with boson field tp and tunneling 
to the leads and back. The quantum kinetic equation (f2"5")) 
is derived for any values of g r and gi; the r.h.s. of Eq. (|25[) 
can be written as the series in powers of g due to the pres- 
ence of Im2?^(r) and B u . The boson distribution B u is 
determined by electron distribution function F d (r) and 
should be found from the solution of the AES-action (fTU| . 

At g 3> 1 the kernel (ImX^) of the collision integral 
of the quantum kinetic equation (|25[) resembles the ker- 
nel of the collision integral in the quantum kinetic equa- 
tion for disordered electron liquid 4 ^— for energy trans- 
fers ui » gS as it is expected^ At g <C 1 the quantum 
kinetic equation (f2~5j) which takes into account the renor- 
malization effects via Im generalizes the kinetic equa- 
tion derived in Ref. [26j in the framework of the orthodox 
theory 4 ^ for sequential tunneling and inelastic cotunnel- 
ing approximations. 



to be equal to the electron-hole distribution function B u 
introduced in Eq. ([2(7|l : 



BUr) = 



Kir) 



2i Imn«(i 



E 9a Jde l-F e d (T)F e a _Jr) 



E 9c Us F?(t)-F?_ u (t) 

y.—l,r 



(28) 



The heat current conservation I? + H = determines 



the equilibrium temperature of the island: 
,( oq ) _ giTi + g r T r 



T 



91 + 9r 



(29) 



A straightforward computation of charge and heat cur- 
rents gives 



gig,- 



4tt7 2^sinh 2 ^ 



(3uj 

/3ui 2 
2 



G8u 




(30) 



Introducing the quantities g',g T , and k' as 
„2 



M K 



we obtain 



9l9r 



h {gi 



9 



ik' 



(31) 



IV. TRANSPORT COEFFICIENTS 

Using quantum kinetic equation (|25|) we are able to 
derive general formulae for all linear response coefficients 
of the SET for any value of g. Voltage (AV = V r - Vfi 
and temperature (AT = T r — Tj) differences across the 
SET cause charge (1^) and heat (1^) currents. Electric 
and thermoelectric transport coefficients are defined as 



_ (G v G T \ (AV 
/(«) J ~ \M K ) \AT 



(26) 



Here coefficients M and Gt (the response of a heat cur- 
rent to voltage difference and the response of electric cur- 
rent to temperature difference, respectively) are related 
via Onsager relation M = GtTA 7 * The thermal conduc- 
tance k is usually defined as k — K — GyTS 2 where 
S = Gt/Gv stands for the thermopower. The electric 
and heat currents in the a-th reservoir can be found as 



2^ 



(27) 



x{(F e d +u (r) Fe(r))BM) 1 + K (r)F £ d +u (r)} . 

The current conservation corresponds to the condition 
If + I£ = 0. It fixes the boson distribution function B^ 



dco Im£>^ 
4^sinh 2 ^ 



flu 
3 



(32) 



We stress that Eq. (|3"2"j) is valid for any value of tunnel- 
ing conductance g. It generalizes expressions for trans- 
port coefficients obtained in Refs. [20ll48l | for g <C 1 to 
arbitrary values of g. 

It is worthwhile to express Eq. (|3"2"j) in terms of the 
tunneling density of states of electrons in the island (see 
Appendix ICj) : 



Vd{e) 



/ Im Vf: \ coth tanh \ — . 

J u \ 2T 2T ) 2tt 



(33) 



Substituting expression (|3"3"|) for the tunneling density 
of states and performing standard integrals with Fermi 
and Bosc distribution functions one can check that the 
results Q32p are allowed to be exactly rewritten in the 
form 




(34) 



Eq. (|34[) for the transport coefficients resembles the cor- 
responding expression in the Fermi- liquid i 47 ' 49 However, 
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contrary to Fermi-liquid, the tunneling density of states 
Vd(e) has strong dependence on energy for e — > 0. In 

general, Vd(e) = ^ vcn (e) + v^ dA {e) where ^ vcn ^° dd (£) is 
even/odd function of e. It can be shown that ^ ven ^ odd (e) 
is even/odd function of the external charge q. Therefore, 
g' and k' are even functions of q whereas g' T is an odd 
function of the external charge. 

For macroscopic samples of ordinary metals, the 
Wiedemann-Franz law provides a universal relation be- 
tween the conductance and thermal conductance. It 
states that the Lorenz ratio C = k/(GvT) , is a constant 
given by the Lorenz number 7r 2 /3e 2 . As follows from 
Eq. (|34p . one can expect the violation of Wiedemann- 
Franz law in the presence of strong dependence of ^(e) 
on electron energy. 

In the case j > 1 one is able to perform perturbative 
expansion in 1 / g and take into account non-perturbative 
corrections. The function lmT> R acquires the following 
form in the equilibrium^ 



lmV R , = - 



T 



2 gE c e< 
l --g Xn ^T 



9 2 E C 

7T 2 T 



-fl/2 



cos 2nq 



1 



2tt 
guj 

g 2 E c _ 



g 3 E c _ 



s/2 



cos 2irq 



9/2 



cos 2nq- 



g 2 E c 

T 



irT 

- g/2 sin27rg ( 5{u 



4ir 2 T 2 
2T 



4tt 2 T 5 



(35) 



Here, function lo8(uj) can be understood as Ima/[7r(w + 
a + iO)] where the limit a — > should be performed at the 
very end of all calculations, (this calculation can be, e.g. 
integration over a;). The non-perturbative in 1/g cor- 
rections (exponential terms exp(— g/2)) come from Kor- 
shunov instantons^i of the AES-action. Then by using 
Eq. we find from Eq. ([MJ) 



k' 



g =g 



9t 



2 In 



gE c g 3 E c 



-9/2 



T 6T 

2g*E /a / it* 

irT V 12 
4 2g 3 E c 

3 7T 2 T 



cos 2irq, 



cos 27rq 



sin 27rg, 

TT 2 _ 

3 



(36) 
(37) 
(38) 



The result for g' has been obtained in Ref. [37|. Equa- 
tions (|37|) - (j3"5)) are new and valid for temperatures T 
g 2 E c exp(—g/2). We emphasize that g' T has only non- 
perturbative instanton) contribution. The same holds 
for the thermopower: 

^-^(l-£)cxp(-f)sm2^. (39) 

At g 1 the violation of the Wiedemann-Franz law is 
weak and the Lorentz number is given as 



£ = — 

3e 2 



4 2/ 9 
1+ % + 3( 1 "^ 



9 2 Ec c/2 
T 



cos 2irq 



Due to the presence of the non-perturbative contribution, 
the Lorentz number is temperature dependent and oscil- 
lates as a function of the external charge q. Eqs. (|39|) 
and P0)l constitute one of the main results of the present 
paper. 

In the strong coupling regime g 1, Eq. (|32[) supple- 
mented by the proper expression for lmD R (cf. Eqs. (|75[) 
and (|114[) ) results in exactly the same expressions for the 
transport coefficients as obtained in Refs. [l7lfl9ll20ll48j . 
We refer a reader to these works for details. 



V. RELAXATION OF ELECTRONS IN THE 
ISLAND, WEAK COUPLING REGIME g > 1 

Next, we want to illustrate the ability of quantum ki- 
netic equation ([2"5]) combined with fine field-theoretical 
scaling of essential physical quantities. We consider the 
problem of relaxation of electrons in the island towards 
the equilibrium due to the tunneling to the reservoirs and 
back. There are two possible scenarios. The first one can 
be refered to as a quasi-equilibrium regime. The elec- 
tron distribution inside the island is given by the Fermi- 
function but with non-equilibrium temperature which 
slowly relaxes to its equilibrium value. The second sce- 
nario is fully non-equilibrium regime when electron dis- 
tribution is arbitrary. Which scenario persists depends 
on the ratio te / T ee where te stands for the energy relax- 
ation time due to tunneling mechanism and r ee for the 
energy relaxation time due to electron-electron interac- 
tion in the island. The non-equilibrium regime persists 
provided te -C T ee and the quasi-equilibrium regime is 



(40) 



possible if te 3> r ee . We will argue below (see Sec. IVIip 
that both scenarios are possible. 

There is one more relaxation time involved: trc which 
determines relaxation of the electric charge on the island. 
In the weak Coulomb blockade regime, trc is given by 
the following clasical estimate: trc — 2tt /gE c . As we 
shall see below, te S> trc- Therefore, it is allowed to 
assume that at first there is quick relaxation of the elec- 
tric charge on the island and, then, slow relaxation of 
the electron distribution function or temperature towards 
the equilibrium. Technically, it means that initial elec- 
tron distribution function F?(0) satisfies the constraint 
fde[F e d (0)-Fn = 0- 

As was discussed in the Introduction, the renormaliza- 
tion of physical observables drastically changes the relax- 
ation dynamics of the system. Therefore, before solving 
kinetic equation we need to establish the scaling of a 
theory's coupling constants under non-equilibrium con- 
ditions. 



A. Renormalization of AES-action at g ^> 1 

The AES-action is renormahzed due to its nonlinear 
form. In the equilibrium case renormalization of the ac- 
tion is well-known (see e.g., Ref. [l|). In our case non- 
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equilibrium makes the problem non-trivial. As in equi- 
librium, we expect the necessary scaling of the coupling 
constant g. The additional question that inevitably arises 
is whether the structure of the kernel of AES action (the 
components of polarization operator n w in Keldysh space 
in Eq. (|16|l ) is changed due to renormalization? The 
details of the calculation are presented in Appendix [D] 
We prove that the structure of the bare action is fully 
restored, the kernel of the AES action being intact dur- 
ing renormalization group (RG) procedure. The coupling 
constant rcnormalizes according to 



9(A) = 5(A) 



duj. 
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(41) 



Here high energy scale A is naturally set by the first term 
in Eq. (|14[) : A ~ gE c . To demonstrate that the integral 
in Eq. (|41[) is indeed logarithmic we explore the behavior 
of the integrand at lo — > oo. It is straightforward to get 
the following asymptotic for function at lo — > oo: 



B u = sgnw + 8Bu, 
SB U = OFj-sgn^) 



(42) 



F«)d E F?de. 



B. Non-equilibrium regime 

The relaxation problem is formulated as follows. At 
t — the island is heated and some electron distribution 
function Ff (0) is created. The characteristic energy Ed of 
electrons in the island is larger than temperatures of the 
reservoirs, which are kept fixed and equal to each other 
£d > T r = T[. The system is released and the island 
is cooling down due to the tunneling of electrons to the 
reservoirs and back. 

Performing expansion to the second order in boson 
fields <p c ,q, one straightforwardly finds (see Appendix |D|) 



r> 2it5(lo) 
ImV* = — 



1 - 



1 



^dLO) - *L. 

lo J gui 



(47) 



We mention that this result generalizes the pertur- 
bative (independent of q) part of Eq. (f3"5j) to the non- 
equilibrium case. With the help of (|47[) one can compute 
the collision integral in the r.h.s. of the quantum kinetic 
equation (|25[) and obtain 



d T F e d (r) = -^(2?(r) 



2tt 



G(r) = g - 



Mr) 



(48) 



duo. 



We expect that any physical distribution function obeys 
the condition F^ — > sgn e at e — > oo. Then 



lim 8B U 



0. 



(43) 



Therefore, the high-energy asymptotic of function B u is 
given by sgne as in the equilibrium. This way, the loga- 
rithmic behavior of integral in Eq. (|41j) is ascertained. 

To get the renormalized action one has to integrate out 
all the frequencies down to the lowest scale ljq, at which 
the RG stops. This energy scale can be determined as 



8B(uo) ~ 1. 



(44) 



Let Ed be a characteristic energy scale of the island distri- 
bution function (the scale at which electron distribution 
function F^ becomes almost equal to sgne). Then one 
can easily check that the following estimate holds (see 
Appendix [Dl for elaborate details) 



w - max{e d , T r , T t }, 



(45) 



where T r ,T; are temperatures of the reservoirs. Energy 
scale loq serves as a natural lower cut-off, A = ujq, in the 
RG procedure. Finally, we find 



g(uo) 



2 In 



gE c 

L0 



(46) 



In the equilibrium, Ed = T r = T\ = T and one finds 
loq = T. Eqs. (|45 |l -(|46 p describe renormalization of the 
AES-action under non-equilibrium conditions. 



Here we neglect last term in Eq. (|4T|) for the follow- 
ing reasons. It gives contribution of the order of unity 
whereas the first term in Eq. (|47f involves J dco B u /ui ~ 
In gE c /£d 7S> 1. Although, Eq. (|4"5)l has a quasi-elastic 
form, in fact, it is highly non- linear equation: G(r) in- 
volves information about the electron distribution at all 
energies. 

As was shown above, the quantity Q (r) has meaning of 
the renormalized coupling constant of the theory. Sim- 
ple algebra leads us to the differential equation for the 
function Q(t)£1 



2tt 

. du oj 

Qr =9-1 — COth— . 



(49) 



The solution reads 
G(r) 



g(o)g r 



g(P) + (g r -g(P))e- 



Sri, 



(50) 



Now by using this result we integrate Eq. (|4"5]l and obtain 
the evolution of the electron distribution function F?: 



F*{t) = FJ+ (F e d (0)-F e r 



g r s 

exp , — r 



- 1 



(51) 



Eq. (|51[) demonstrates energetically uniform relaxation 
of the electron distribution function. This fact is a direct 
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consequence of the quasi-elastic form of the kinetic equa- 
tion (|48[l. However, due to renormalization effects the 
form of the relaxation law is different from the exponen- 
tial one. 

Let us define the characteristic energy e d as e 2 d — T ? 2 = 
(3/tt 2 ) J dee(F^ - F 6 d ) such that e d = T d in the quasi- 
equilibrium case = tanh(e/2T d ). Then, in the case 
£d(0) 3> T r and at not too long times t -C 2ii/8G ri 
one finds from Eq. (|?T|) that the characteristic energy 
decreases according to the power-law: 



■ >E ch (n,q) 



£d{r) = e d (0) 



SG(0) 

2ir 



-1/2 



(52) 



C. Quasi-equilibrium regime 



In the quasi-equilibrium regime, 



we 





need to take 

into account the collision integral /| ee ' ) due to electron- 
electron interaction in the islandi 43 ' 44 As this term is 
added to the r.h.s. of Eq. (|2"5]) . it makes the electron 
distribution to be the Fermi-function. Multiplying 
both parts of Eq. (|2~5"j) by e and integrating them over 
energy, we obtain the following equation (the well-known 

r (ee) 



identity J deelg = is used): 



dT 



g<[ (T 2 

2n {Id 



(53) 



Here we use the leading (classical) part of Eq. (|4"T)) 
(lmT>% = -2ttS(lj)/B u ). Equation (|53J| yields standard 
exponential relaxation towards the equilibrium. In the 
limit Td » T r (£(0) ^S> Q r ) it is also possible to compute 
collision integral using the entire one-loop expression (|4"T)) 
of its kernel. Naturally, one-loop correction reveals itself 
in the logarithmic renormalization of g in the r.h.s. of 
Eq. ([53]). By using Eq. (|47j) , we find 



dT 



9(t)6 
2ir 



G(r)=g-2\n 



gE c 

cT d {r) 



(54) 



where c is a numerical constant of the order of unity 
which docs not influence final results. The solution 
of ([54]) reads: 



Td(r) =T d (0)exp 



2tt 0(0) 



G(0) 



(55) 



The condition in the second line of Eq. (|55|) implies that 
solution holds for not too long times at which T d (r) >• T r 
(G(t) 3> Gr)- The logarithmic renormalization of the 
conductance changes the character of temperature relax- 
ation. At long times 2ir/6 < r < (2k /S) ln(?(0)/</ r , the 
cooling of the island slows down in comparison to the 




FIG. 3: (Color online) Charging energy E c h = E c (n — q) 2 as 
a function of gate charge q. 



standard exponential decay which is developed at short 
times t <§C 2ir/5: 



T d (r) = T d (0)e~ — 



(56) 



It is instructive to compare the relaxation of temper- 
ature in the quasi-equilibrium regime and the character- 
istic energy e d in the non-equilibrium regime given by 
Eqs. (|56"|) and (f52"|) at times r -c 2ir/8Q r , respectively. 
While the former demonstrates exponential behavior, the 
latter decreases in accordance with the power-law. 



VI. RELAXATION OF ELECTRONS IN THE 
ISLAND, STRONG COUPLING REGIME, g < 1 

In the strong coupling regime there are two possible 
scenarios for relaxation of electrons in the island of the 
SET. The first one persists if te <C r ee . In this non- 
equilibrium case the carriers inside the island do not have 
time to thermalize and to form the Fermi-distribution 
with some temperature. In this case the the time evo- 
lution of distribution function itself becomes the main 
objective. This task is solved in section IVI CI below. The 
second scenario develops in the opposite limit, te 3> T ee . 
Namely, the relaxation rate due to electron-electron in- 
teraction inside the island is much faster than the rate 
due to electron tunneling through the contacts. Thus, 
the temperature of carriers in the island becomes a well 
defined characteristic of a system. Consequently, the re- 
laxation of the island's temperature will be the focus of 
our analysis in section IVI Dl 

As in the previous section we shall assume that the 
electric charge on the island quickly relaxes and only then 
slow relaxation of the electron distribution or tempera- 
ture starts. In the strong Coulomb blockade regime this 
picture is well justified since 1/trc — gmax{T, A} 3> 
1/te. 

We concentrate on the most interesting case: the vicin- 
ity of a degeneracy point: q — k + 1/2 where k is an in- 
teger. Following Ref. [53| . the hamiltonian (dl"© can 
be simplified by truncating the Hilbcrt space of elec- 
trons on the island to two charging states: with Q = k 
and Q = k + 1 (see Fig. [3]). The projected hamilto- 
nian then takes a form of 2 x 2 matrix acting in the 
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space of these two charging states. Denoting the devi- 
ation of the external charge from the degeneracy point 
by A: q = k + 1/2 — A/(2E C ) we write the projected 
hamiltonian as^ 

H = H + H t + AS Z + A_ + El (57) 



where Hq is given by Eq. ^ and 

H t = t ka a\d a S~ + h.c. 



(58) 



Here S z , = S x ± iS' y arc ordinary (iso)spin 1/2 op- 
erators. 



A. Non-equilibrium pseudo-fermions 

To deal with spin operators it is standard to use 
Abrikosov's pseudo-fermion technique.— We introduce 
two-component pseudo-fermion operators ijj^, ip a such 
that 



(59) 



The out-of-equilibrium pseudo-fermions were tackled be- 
forei 55 i 56 As usual, one introduces Keldysh contour, dou- 
bling the number of fermions. The system is out of equi- 
librium and one has to be very cautious. The distri- 
bution function T e of pseudo-fermions is not known a 
priori. Rather, it is to be defined self-consistently from 
corresponding kinetic equation. Pseudo-fermions are also 
subject to constraint on their number: 



i 



(60) 



Thus the state of a system ought to be projected on the 
state with M = 1 at any instant of time. The operator of 
particle number is conserved by Hamiltonian (|57[) - (|58|) . 
Consequently, the operator of projection on to physi- 
cal subspace Af = 1 commutes with hamiltonian too. 
It means that the projection on to physical subspace is 
needed at a single point of Keldysh contour only. We in- 
sert the factor exp(r) J2 a V'aV'a) into density matrix and 
take the limit rj — > — oo at the end of any diagrammatic 
calculation. Then 



(0) - Km 

ri-t-ao {Al)p{ 



(61) 



Provided the operator O has zero expectation value in the 
sector with zero pseudo-fermion number, M = 0, Eq. (|61[) 
can be simplified as 



(O) = lim 



(0) P i 



r/^-oo (A/") 



pf 



(62) 



The dissipative action is to be rewritten in the Keldysh 
representation. We plug representation (|59[) into the 



hamiltonian (|57[) and integrate out electrons in the lead 
and the island. This leads to the following effective action 



5' 



dttp (idt 



(J, A 



r))ip 
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^(t)7i<r_V'(t)n iJ -(t,t')^(07j^+V'(*')d*d*' 



(63) 



Here a stand for Pauli matrices, a± = (a x ± ia y )/2, and 
7i = 7 9 = ( ? n ) , 72 - 7c = (I ? ) (64) 



are matrices in Keldysh space. The pseudo-fermion op- 
erators ?p are understood as vectors in the tensor product 
of isospin and Keldysh space, stands for the matrix 
of polarization operator (|TB]) - (|T7J) . Next, we write the 
Wigncr transform of the quantum kinetic equation for 
pseudo-fermion distribution function 

-id T F e (r) = Sf (r) - Ef (t)T £ (t) + ^(r)-F e (r). 

(65) 

Here as before, we neglected all derivatives with respect 
to slow time r. All functions entering Eq. (|65|) are un- 
derstood as matrices acting in the isospin space. From 
the appearance of Eq. (|6"5)) we conclude that character- 
istic relaxation time of pseudo-fermion distribution func- 
tion F s (t) is T p f ~ l/(<7max{A, T}) and is much shorter 
than te- 

It allows us to consider pseudo-fermions to be in the 
stationary state. Then the l.h.s. of the kinetic equa- 
tion (|65[) can be omitted and we obtain the equation for 
pseudo-fermion distribution function: 



Mr) 



2iIm£| ? (T) 



(66) 



With the help of Eq. (|fJ3"| we write down equations for 
the pseudo-fermion self-energies (Fig. [7]): 



E + (M')^E n ^G 



tt'-li 



(67) 



Here Gtt',a stands for the pseudo-fermion Green func- 
tions corresponding to the first line in Eq. (|6"3"]l and £± 
are matrices in the Keldysh space. We will need the ex- 
plicit expressions for their Wigner transforms: 



ig f duj 



= ~J ^ImnJlmGf^^.^j;-^^ - a}, 
lm^ a = - 9 -J ^Imn>G E fl +TOii _,{B u -<rJ^), 
9 J ^Imn«ReGf +CTW> _ CTJ B w . (68) 



ReS" =-■- 
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Here a stands for ± and 



Ac 



(69) 



Combining Eqs. and wc find the following equa- 
tion for the pseudo-fermion distribution function 



-a(e+^+fj)« 



(70) 



where F a = -?X<t/2-7; ■ P m g§ m S £ = Acr/2 — rj we arrive 
at the closed equation for F a : 



(71) 



Now we need to investigate asymptotic properties of 
functions T a when pseudo-fermion chemical potential 
ij — > — oo. It is natural to expect that the equilibrium 
result 



lim T a = 1 

Tj — V — OO 



(72) 



survives in the non-equilibrium. As one can check this 
assumption satisfies Eq. (fTTj) . 

In order to solve the quantum kinetic equation (|25|l. 
we need to compute Iml?^ in the strong coupling limit 
g < 1. With the help of Eqs. (J3HJ) and (JHHJ), one easily 
finds in the zcroth order in g: 



li\\V, 



R 

aj.pf 



^ImGf + -ImGf + l 



(73) 



Now we express the physical correlation function through 
pseudo fermion one Im2>^ pf . By using the following 
zcroth order in g result for the pseudo-fermion number 

(ao^ =£ / t l ^ G U^ - 1 ) = 1 - 



2tt 

we obtain 

ImD^ = —n5(uj + A) lim 



(74) 



T- + 1 
s-— oo 2B_a + 1 ~~ -F 



tt5{lo + A) 



(75) 



Eq. (|75[) is the generalization of the equilibrium result 
for correlation function Im2>^ (see Rcfs. gilj) 

over 

the non-equilibrium case. 

Next, by using Eqs. (|55j) and (f5U)) , we find in the zeroth 
order in g: 



££pf = "2* / ^ImGf +Wi _ImG« + (l-J-- +a) J?- 



cfe 



ttzci(w + A)(J"+J- - 1) 



(76) 



Expressing the physical correlation function through 
pseudo-fermion one { , we obtain 



2?5 = -27ri(5(a; + A). 



(77) 



This result implies that the boson distribution function 
Bu, is determined by B u in the same as in the weak cou- 
pling regime, 



(78) 



Before proceeding with the solution of the quantum 
kinetic equation we prefer to perform one-loop renormal- 
ization of the theory. This is done to sum up all large 
logarithmic corrections (which otherwise arise in pertur- 
bative analysis) and absorb them into renormalized phys- 
ical constants of the theory. 



B. One- loop structure of the pseudo-fermion 
theory 

In this section, we establish the out-of-cquilibrium gen- 
eralization of the scaling of fundamental parameters in 
the pseudo-fermion theory (the gap A, the coupling con- 
stant g), the Green's function and the average pseudo- 
fermion density (Af) p f. We expect that the action (|63|) 
can be renormalized with only one scaling parameter Z 
like in the equilibrium case [57H59| . This is indeed the 
case and the obtained renormalized structure of the the- 
ory is a natural generalization of the equilibrium one. 
The renormalized pseudo-fermion Green's function be- 
comes 



Z(X) 



e-£, a ±igT a {e)' 



la = -r\ + oA/2, (79) 



where 



A 



f>' (80) 



See Appendix [E] for details of the computation. It is 
straightforward to check that coupling constant g and 
gap A are renormalized according to: 



g = gZ 2 (X), A = AZ 2 (A). 



(81) 



To complete the renormalization picture we need to es- 
tablish the scaling dimension of the pseudo-fermion num- 
ber (Af)pf- In complete analogy with Ref. J57|, (A%/ 
happens to have no renormalization 



(N) pf = (N) P f 



(82) 



For completeness we present the rigorous proof of 
Eq. (|82j) via Callan-Symanzik equation in Appendix [E] 
The integral in Eq. (|5U)) runs over frequencies E c 
^> luq = max{T r! Ed, A}. The energy scale ujq deter- 
mines the natural scale at which the RG procedure has 
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to be stopped. The Green's function (1751) acquires the 
width 

Me) = ^(e - + B e _ f _J, (83) 

where = J 7 / . Therefore, the renormalized physical 
correlation function becomes 



Im^ = -Z 2 (A)^ (CJ + A) 



-2niZ 2 (\)5(uj + A) 



(84) 



C. Electron distribution relaxation in the island 

In this section we consider the relaxation in the non- 
equilibrium case, te <C T ee . We focus on the most in- 
teresting case of the Coulomb peak: A = 0. Then the 
quantum kinetic equation (|25[) is greatly simplified (cf. 
Eq. gHD): 



Q(r) 



1 



2tt 2 



(85) 
(86) 



Here, we stress that the kinetic equation ([55")) - ([81)1) is of 
the infinite order in the electron distribution function 
on the island. Indeed, A involves Ff via electron-hole 
distribution function B u . 
The formal solution reads 



F e d ( T )=F^ + (F e d (0)-F^e^> 



(87) 



The function Q (r) obeys the differential equation 



d T Q(r) = -^5 2 (r) 



5(r) 



is given as 



/ Qr 



f 



Qr 



The solution of Eq. 

sg r r _ 

2tt _/ \G(0)J J \G(t) 
f(z) = z + ln(l-z) 

By using the relation 

SQ r T Q r 



(89) 



r[ T dT' SQ r 

6 Jo ^ G{T) = —n 
which follows from Eq 

F 6 d (r)=^; + (F £ d (0)-F:)exp 



Qr 



Q(0) Q{t) 
we obtain 

SQ r T Q r 



(90) 



2tt 



Q(0) Q{t) 
(91) 



Since Eq. f|89[) can not be solved analytically with re- 
spect to Q(t) it is instructive to investigate limiting cases. 



Let us assume that the effective energy of electrons in 
the island Ed 3> T r such that Q(0) Q r - Then, expand- 
ing f(z) in the series in z, we find 



F^r)=F-+{F e d (0)-F: 
Q r [ 



x exp 



SQ r T 



Q(o) V £ 2 (°) 



(92) 



Eq. ([9"2"[) is valid provided Q(t) 3> Q r , i-e., for not too 
long times: r <C it / 5Q r . It is worthwhile to mention that 
standard exponential relaxation 



F £ d (r) = Fl + (F*(0) - FJ) exp 



6G(0)t 
2tt 



(93) 



occurring at short time r <C TtQ r / (SQ 2 (0)) transforms 
into regime of slower relaxation at intermediate time 
^a,7(^ 2 (0))<r< w/(5Q r ): 



F«{t) = Fl + (i^(0) - Fj) exp (- 



SQrT 



(94) 



At longer time t ^ 7r/(<5^ r ), function Q(t) becomes al- 
most equal to Q r and we find again the regime of standard 
exponential relaxation: 



F e d (r)=F; + (^(0)-F;)exp 



SQ r T 



2tt 



(95) 



The same exponential relaxation as given by Eq. (|95[) 
holds if the effective energy of electrons in the island £</ 
is slightly larger than T r such that Q(0) — Q r -C Q(0), Qr- 



D. Temperature relaxation in the island 

Now we investigate the relaxation in the quasi- 
equilibrium case, te ^> T ee . 



1. Coulomb peak, A = 

We start from the regime of the Coulomb peak: A = 
0. In the quasi-equilibrium regime, one needs to add 
to the r.h.s. of Eq. (|85j) the collision integral /i ee ' due 
to electron-electron interaction in the island. It is this 
term that makes the electron distribution to be a Fermi- 
function. By using the well-known identity J deel^f^ = 
0, we obtain the following equation: 



dT 



Q(t)5 
2tt 



T 2 d {T) - t; 



(96) 



where Q(t) is given by Eq. (|86[) . In the quasi-equilibrium 
case, we can not derive closed equation for Q(t) as it was 
done in the non-equilibrium case due to the presence of 
additional term /J ee ' ) in the r.h.s. of Eq. ([55]) . 
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Assuming that T<j(0) >• T r we can estimate A with 
logarithmic accuracy as A = In E c /T d . Then, we find 
from Eq. (J96J) 



Q(t) = G(0) 



SG 2 (0)t 



2tt 3 



-1/2 



and 

T d {r) =T d (0)cxp 



5(0) 



<5£ 2 (0)r 



2tt 3 



(97) 



(98) 



The solution (f9"5)l is valid provided the condition Trf(r) ^> 
T r holds. If T d (0) > T r exp(7r 2 /£(0)), then the exponen- 
tial relaxation 



T d (r) = T d (0)exp 



(5g(0)r 

47T 



(99) 



developing during initial period r <§C 27r 3 /[(5^ 2 (0)] trans- 
forms into regime of slower relaxation at intermediate 
time: 



T d (r) = T d (0)exp 



7r5r 



(100) 



2tt 3 2, 2 T d (0) 

We mention that in this regime the temperature re- 
laxation is independent of the quantity G(0) which de- 
termines the SET conductance. In the opposite case, 
T d (0) < T r exp(7r 2 /£(0)) the temperature T d {r) evolves 
according to Eq. flM]) for r < {4m/8Q{G)) lnT d (0)/T r . 

At longer times r > (47r/<5£/ r ) \nT d (0)/T r the temper- 
ature T d {r) becomes of the order of T r : T d (r) — T r <C T r 
and we find the standard exponential relaxation: 



T d (r) = T r + (T d (0) - T r ) exp - 



SQ r r 

47T 



(101) 



Evolution of the temperature of electrons in the island 
is presented in Fig. We mention universality of the 
relaxation at long time when the difference between the 
electron distribution in the island and in the reservoirs 
becomes small. In non-equilibrium te -C r ee and quasi- 
cquilibrium te 3> r ee regimes the relaxation is exponen- 
tial with a rate of the order of SQ r . The same exponen- 
tial relaxation as in Eq. (jlOip holds if the temperature 
of electrons in the island T d (0) is slightly larger than T r , 

T d (o)-T r <r d (o),T r . 

It is worthwhile to mention that there is a parametric 
region of time domain g r /(Sg(Q) 2 ) < r < l/(5£(0) 2 ), 
when the relaxation of the distribution function in 
the non-equilibrium regime is much slower ln(i^(r) — 

F I)/( F e(°) ~ F e) v 7 ^ th an the relaxation of the 

(Fermi) distribution function in the quasi-equilibrium 
regime, i.e., relaxation of temperature, ln(T d /T d (0)) ~ 

— T . 



T d (t) - T r 




e 2 (o) 



In' 



T d (0) 



St 



FIG. 4: The dynamics of temperature relaxation, j < 1, 
Td(0) > r r exp(7r 2 /e;(0)). Here, &(t) ~ -Q{Q)5t, &(t) ~ 

-VSi, ~ -GrSt. 



2. Coulomb valley, A > T d (0) 

Now we consider the relaxation of the electron tem- 
perature on the island in the regime of Coulomb valley, 
A ^> T d (0). By using Eq. ((84)) . we rewrite the quantum 
kinetic equation (|2"5]l as 



d T F a e 



4tt [ 



1 - F*F 



e+A 



B-A 



(102) 



We remind that we consider the quasi-equilibrium 
regime. Then we need to add to the r.h.s. of Eq. (| 102[) 
the collision integral /I 66 ' which describes scattering due 
to electron-electron interaction in the island. In what 
follows we assume that the condition T d (0) ^> T r holds. 
With the help of the following results 



dse(l-F e d F^J=T^gnco 



4|w| 



-41i 2 (-e-l w l/ T ") + pin 1 + e-^ T * 
T d V 

B u = — - In | 2 cosh — 



(103) 
(104) 



which are valid for T r < T d (li 2 (z) = E^Li^A* 2 
denotes the polylogarithmic function), we obtain from 

Eq. (una) 



dr 47T 3 



A 



(105) 
(106) 



We can estimate parameter A with logarithmic accuracy 
and find A = lni? c /A since the temperature of electrons 
in the island T d <C A. Therefore, both g and A are 
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independent of r. Integration of Eq. (|102[) yields 

h(z) = e z /z-Ei(z). (108) 

Here Ei(z) = — dtexp(—i)/t stands for the integral 
exponential. By using the asymptotic h(z) = exp(z)/z 2 
at z ^ 1, we obtain 



q{t) = g(o) 



In 1 



T d (0) 
A 

3(5Ag(0)r 
+ 47r3T d (0) 



3<5Ag(0)r - 
4^(0) , 
l 



(109) 



and 



T d (r) = T d (0) 



A 



In 1 



36A£(0)t 
47r3T d (0) 



n -1 



(110) 



The results (|109[) and (|110[) are valid at not too long 
times 



4^ 3 T d (0) / A 

T<< 3^(bT cxp U 



(111) 



As expected, due to the exponentially small SET conduc- 
tance in the sequential tunneling regime, the temperature 
relaxation is very slow, namely, logarithmical. Therefore, 
it is instructive to consider contribution to the tempera- 
ture relaxation due to the electron co-tunneling. 



3. Inelastic cotunneling regime 

As known very well, due to exponential suppression of 
the sequential tunneling mechanism deep in the Coulomb 
valley, Td <C A, the higher order process of inelastic 
cotunneling dominates the transport.— Contrary to the 
case of sequential tunneling, the cotunneling contribution 
to the collision integral in the r.h.s. of Eq. (|2"5)) comes 
from frequencies of order u ~ Td -C A. 

In the pseudo-fermion technique the inelastic cotun- 
neling is revealed as the broadening of delta-peaks in 
the imaginary part of the retarded and advanced pseudo- 
fermion Green functions After taking into account 
Eq. (|83|) , the integrand in Eq. ((73)) becomes of a complex 
pole structure. There are two pairs of proximal poles 



C+±zgr + (e + ), 

S--w±igT _(£_). 



(112) 



There is an additional series of Matsubara-type poles re- 
sulting from distribution functions J-" e + and F~ +La . They 
lead to logarithmically divergent sums. The latter are 



controlled by the renormalization scheme. In our case all 
leading logarithms are absent. They have already been 
absorbed into renormalizcd constants g and A by the 
proper choice of reference energy scale. Thus we can 
omit all divergent sums over Matsubara frequencies. Ex- 
panding in the w/ A we obtain 



gZ 2 uj T + + J- 
A 1 



, M«|A| 



(113) 



Next we use the same arguments that led us to leading 
order espression ([75)1 . The function Im T>^ then reads 



Im£> 



R 



gz 2 uj 

4tt A 2 ' 



M « |A| 



(114) 



Using Eq. (|114[) . we rewrite the quantum kinetic equa- 
tion (1231) as 



d T F e d 



g 2 S 



16tt 3 A 2 



dui £ 



(F:_ u -F e d )B u + l- F*FI_ 



(115) 



We remind that we consider the quasi-equilibrium 
regime. We mention that Eq. (|115[) coincides with the 
kinetic equation derived for the cotunneling regime in 
Ref. HI]. Then we need to add to the r.h.s. of Eq. (fTT5)) 

the collision integral /| ee ' ) which describes scattering due 
to electron-electron interaction in the island. 
In the case of Td — T r < T r , we obtain 



dTd 
dr 



35 Q r 

5-7T 



(T d - T r 



(116) 



where Q r = g 2 T 2 /(6A 2 ) stands for the equilibirum SET 
conductance in the cotunneling approximation. In the 
opposite case Td S> T r , by using Eq. (|115p . we find the 
following equations: 



dTd 
dr 



S(r) 



20?r 
9 2 T 2 d (r) 
6A 2 



Ur), 



(117) 
(118) 



which govern the temperature relaxation. It is worth- 
while to mention that if one substitutes Q r by Q(t) in 
Eq. (fTT6)) then it becomes similar to Eq.(4) of Ref. [38| 
for V — and in the absence of phonons. However, due 
to different numerical coefficients in the right hand side 
of Eqs. (|116|) and (|117|) such substitution is impossible 
even on the level of interpolating expression. Therefore, 
in the case Td — T r ~ T r one needs to solve Eq. (|115|) 
numerically. 

Though, Eq. ()116p leads to the standard exponential 



relaxation, Eq. (|117[) results in the relaxation according 
to the power-law. Indeed, solving Eqs. (|117l) - (|118|) . we 
obtain 



q{t) = g(o) i 



3(5g(0)r 
IOtt 



(119) 
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and 



T d (r) = T d (0) 



3ag(o)i 

IOtt 



-1/2 



Eqs. (|119p and (|120[) are valid at times 

IOtt T 2 (0) 



T < 



3(J0(O) T r 2 



(120) 



(121) 



VII. DISCUSSIONS AND CONCLUSIONS 

We have studied the relaxation dynamics of the SET 
under essentially non-equilibrium conditions. The lan- 
guage of kinetic equations happened to be the most 
adequate for this task. Analytical results are pro- 
cured in the limiting cases of weak j » 1 and strong 
g < 1 Coulomb blockade. All relaxation equations (see 
Eqs. ([48 ]l .([85 )l .(|96 |l .lfT05 |l ) obtained in the course reveal a 
pleasant generality. Namely, 



x d ~-sg(x d )(x d -Xr). 



(122) 



Here X d is a relaxing physical quantity (temperature, dis- 
tribution function), and Q(X d ) is conductance of a SET 
which depends on X d . Equation (|122ll has a transparent 
intuitive interpretation. Namely, the characteristic time 
scale determined by the r.h.s. of Eq. (|122p is simply a 
dwell time of a particle inside the metallic island^ i.e. 
t^ 1 ~ Q5. The inverse dwell time can be also estimated 
as the ratio of the thermal conductance k to the heat 
capacitance of the island. The latter is proportional to 
T d /S. The generality of Eq. (|122[) is, however, deceptive 
as it leads to drastically different evolution of physical 
quantities over time in the case of small and large values 
of g. 

In the course of all our analysis we generally discarded 
the influence of electron-phonon (e-ph) interaction. The 
reasoning behind this is as follows. The e-ph scattering 
rate was well studied for 2-dimensional electron gas with 
disorder^ The following estimate has been found 



'c-ph 



8.3 x 10 s T 3 [s^K 



(123) 



The electron-electron (e-e) scattering rate in mesoscopic 
systems is widely studied as well (see, e.g.[62|). For small 
diffusive electron systems and for T <C -&rh there are two 
parametrically different situation a 62 ' 63 



T 2 6 _ . 

^Th 



(124) 
(125) 



where £d = {kpl) T=s /kF and L stands for the size of 
the island. Equation (|125[) is a typical Fcrmi-liquid ex- 
pression coming from large momenta of the order of the 



I 1 1 

quasi-cquilibriuni / I 

^ non-equilibrium 

/ ; 



'=©" 



FIG. 5: Schematic diagram of different regimes for j > 1. 
The non(quasi)-equilibrium regime dominates in (un)shaded 
region. 



max{T, A} 




FIG. 6: Schematic diagram of different regimes for g <C 1. 
The non(quasi)-equilibrium regime dominates in (un)shaded 
region. At A = 0, 7 = ^fg and 7 = g in the cotunneling case 
for T < A. 



inverse screening length. The upper one comes from mo- 
menta k ~ 1/L and of diffusive origin. Let us address 
the question which kind of dissipation dominates in var- 
ious parametric regimes. The relaxation due to electron 
tunneling can be roughly estimated as 



1/te ~ Q5. 



(126) 



By comparing Eqs. (|124[) . (|125[) and (|126[) . one can see 
that both quasi-equilibrium and non-equilibrium regimes 
can occur for j > 1 and j « 1. The non-equilibrium 
regime prevails for j » 1 while the quasi-equilibrium 
one dominates for g <C 1 (see Figs. [5] and [6]). 

To estimate e-e scattering rate we use some experi- 
mental data taken from the experiment by Pasquer et 
aim where they studied the Coulomb blockade effects 
in a small island of two-dimensional electron gas. The 
experimental data were as follows: the level spacing 
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S w 85 mK, the Fermi energy Ep ~ 47 K, the elas- 
tic mean free path I « 15 fim, the size of an island 
L w 1 /im. This allows us to estimate the Thouless 
energy £/rh ~ 8 if and e-e relaxation rate 

r" 1 « 1.7 x 10 8 T 2 [if^s -1 ]. (127) 

The typical temperature of the contemporary mesoscopic 
experiment is T < 100 roif. As we see, with lowering 
temperature the e-ph scattering rate decays faster then 
the corresponding electron-electron (e-e) rate. On the 
other hand for the same metallic island the typical relax- 
ation rate due to electron escape to reservoirs is 

r^ 1 ~ g x 10 s [s- 1 ] (128) 

Estimates (|123[) - (|128[) show that for all relevant exper- 
imental temperatures the phonons are frozen and e-ph 
interaction can safely be discarded. Next, comparing 
estimates (|127p and (|128[) we conclude that varying g 
two different parametric regimes explored in this paper 
can indeed be realized in the experiment. Namely, fully 
non-equilibrium regime is realized when g is large enough 
and electron distribution function is arbitrary inside the 
island. The quasi-equilibrium regime persists in the op- 
posite limit, when g is small enough. 

In addition to the relaxation of the electron distribu- 
tion in the island due to escape of electrons to the reser- 
voirs which we consider in details above there is another 
mechanism of energy relaxation which is due to interac- 
tion Ui r of electrons in the island with electrons in the 
reservoirs. For a sake of simplicity we assume that the 
typical interaction parameter r s ~ l/ikpas) ~ 1 with 
as standing for Bohr radius. In the case L <C Co, the 
energy relaxation rate due to interaction Ui r of electrons 
in the island with electrons in the reservoirs can be esti- 
mated as 

~ W d U ir (k F )} 2 ^. (129) 

Tee ™ 

Here Ui r (k) denotes the Fourier transform of the inter- 
action Ui r (r). Provided the condition kpd 1 holds 
Uirikp) is strongly suppressed, VdUi r (kp) <C 1, and 

4) « h (130) 

In the opposite case of large island, L 3> Co, and for 
(ij « I the estimate for the energy relaxation rate 
l/fie becomes 

7>>~ U 2 (k L )-Ul(k L )u d U(k L )E^ h °- { } 

Here k^ ~ 1/L and U(k) stands for the Fourier trans- 
form of interaction between electrons in the island k = 
component of which leads to the charging term H c in the 
hamiltonian (Q]). As one can see, both cases of r ee Tee 
and Tee Tee^ arc possible for L 3> Co- 



For d <C L where d stands for the typical size of 
the tunneling junction between the island and reservoir 
Eq. (|13ip can be simplified as 

4t - ( i32 ) 

r e ( r } d Th 

For the experiments by Pasquer et al£^ we estimate the 
Bohr radius as ~ 10 nm and assume typical d to be of the 
order of 100?ito. Therefore, we expect that the regime 
in which the main mechanism of the energy relaxation of 
electrons in the island is due to its escape to the reservoirs 
can be realized in a laboratory. 

To summarize, we have explored heat transport and 
relaxation processes in a SET with large number of tun- 
neling channels over a wide range of parameters. In 
the regime of linear response we obtained analytical ex- 
pressions for transport coefficients (conductance, thermal 
conductance and the response of electric current to tem- 
perature difference) in the entire span of values of g. It 
is possible to shape the general relations for linear re- 
sponse coefficients into Fermi-liquid type form. There 
is however an important difference, namely: the tun- 
neling density of states undergoes dramatic renormal- 
ization due to Coulomb interaction. The latter leads 
to violation of Wiedemann-Franz law: in the g 1 
limit the Lorentz ratio C acquires weak periodic depen- 
dence on gate voltage (the precursor of Coulomb block- 
ade). The method of quantum kinetic equation supple- 
mented with non-equilibrium AES action has allowed us 
to treat Coulomb interaction exactly. We have obtained 
the time evolution of electron temperature (in the quasi- 
equilibrium regime) and the distribution function (in the 
non-equilibrium regime) of a SET island due to particle 
escape to reservoir. The corresponding collision integral 
is always non-local in energy due to inelastic nature of 
tunneling processes: the radiation of plasmon ip always 
accompanies the tunneling event. In general, this leads 
to highly complicated integro-differential kinetic equa- 
tions. Surprisingly we have shown that kinetic equations 
can be reduced to simple differential ones in a number of 
wide parametric regimes, namely: j> 1 (weakly block- 
aded SET) and g < 1 (strongly blockaded SET in se- 
quential tunneling approximation with renormalization 
taken into account). This simplification is achieved due 
to the presence of strong scale separation in the problem 
gS <C Td or gS <C Ed- Indeed, the characteristic frequency 
at which the distribution function in the kinetic equation 
changes is uj ~ g5, while the scale at which the renormal- 
ization due to the presence of Coulomb interaction occurs 
is u > Td or lu > Ed- This separation is that allows us at 
first to treat Coulomb interaction and secondly to study 
evolution of the distribution function. 

Still, quantum fluctuations of charge significantly 
change the relaxation laws comparing to simple exponen- 
tial ones which are characteristic of semi-classical physics 
for g S> 1 and of orthodox theory for g -C 1. The regime 
g <C 1, A > T is dominated by cotunncling process. 
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In the latter case the kinetic equation retains its integro- 
diffcrcntial structure and is to be solved numerically else- 
where. Measurements of the predicted relaxation dynam- 
ics are an experimental challenge. 
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Appendix A: Keldysh form of AES action 

For a benefit of a general reader, we outline here the 
details of the derivation of the Keldysh form of AES- 
action from hamiltonian ([IJ-Q. To get rid of unsuitable 
quartic electron-electron interaction term ((4]) we decou- 
ple it via Hubbard-Stratonovich bosonic field <p(t). After 
that the initial electron operators are gauge-transformed 
according to 

->• <4(*)e Mt) , da(t) da(*)e-*W (Al) 

the action of the system becomes gaussian in fermions: 

S = So ~t~ S c + St, 

s o= I £4(^-4 d) K* 

V" a\(ie t - el')a k , 
iu (A2) 

Sc = J P 2 dt + q J pelt, 

S t = - Jj2 (tk a a\d a e l * + H.c) dt. 

Here, for a sake of simplicity, we consider an island con- 
nected to a single reservoir. Superscript d refers to the 
island and r - to the reservoir. The integrals are un- 
derstood as contour ones and 7 is the Keldysh contour. 
Integrating out fermions we obtain the effective action 
for the bosonic field ip: 



Here, matrices G, T have the following structure in the 
reservoir-island space: 



G = 



G 



k,d 







Got,! 



, T = 





4* f 



tkaX 





y - _L ( X ° X <i 
j2\X q X C J> 



(A4) 



where X c ^ q are defined in Eq. (fT5|) . Expanding S e s to the 
second order in T. wc find 



Scs = 



GTGT 



Sr. 



(A5) 



S c ff 



-iti \niG- 1 +f) + S c . 



(A3) 



This expansion is valid in the limit g c h <?C 1 and iV c h 3> 1 • 
Computing all the traces we recover the dissipative part 
of AES-action in form (|16[) with the polarization operator 
n given by the following general expressions: 

n R ' A (t,t>) = Y a T,\ t ^\\ G U t '^) G tdM 
+ Gi; r R (t',t)G« d (t,t')), 

+ Gf,, r (t', t)G% td (t, t') + G* r (t', t)G* d (t, t')) . 

(A6) 

Provided the density of states of electrons on the island 
and in the reservoir are slow varying near the Fermi en- 
ergy, we can perform the summation over a and k with 
the help of Eqs. ©-([5]) and reproduce the kernel of the 
action in form of Eq. (fTT|) . 

Appendix B: Electron's self-energy 

Here, we present the expressions for electron's self- 
energy to substantiate the derivation of the kinetic equa- 
tion in Sec. IIIII As follows from Fig. [2] a the lowest order 
(in 1/iVch) contribution to the electron's self-energy is 
given by 

act' 

k 

a a' 

x E4W (Gt(M')^(M') + [GtM - Gt-(M')] 

k 

x [V R (t,t')-V A (t,t')] 

(Bl) 

In case of constant densities of states in the island and the 
reservoir, they can be simplified with the help of Eqs. (|6])- 
(|8"J) and, then, written in the form of Eq. 
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Appendix C: Tunneling density of states on the 
island 



The tunneling density of states of electrons inside the 
island is defined via corresponding full retarded Green's 
function of original fermionic operators: 

iG R (t,t') = ^(d+d'+e-^^-^ - d+d'_e- i ^+-' p -^ 
+ d-d'+e-^— O - d_<?_e _< ^— v '- ] ) 

(CI) 

Here, operators d± = d±(t), d± = d±(t') are the 
gauge transformed operators of electrons inside the island 
(see Eq. IA11 ) Subscripts ± correspond to upper(lower) 
branch of Keldysh contour. 

Switching to Wigncr transform we obtain 



(C2) 



G%(e) = -J2 / {G2 td (e + u,)B u 1mVl 
+V A F? +u lmGf(s + u)}^-. 



Then, the tunneling density of states of electrons on the 
island becomes 



du> 



Me) = ~Tm Gf(e) = v d j ImD«{^ " j 2 

(C3) 



Eq. (|C3|) gives the tunneling density of states of electrons 
on the island in a non-equilibrium regime with arbitrary 
electron distribution function Fd- In the equilibrium, it 
leads to the result (|33l) . 



Appendix D: Renormalization of AES-action at 

<7>1- 

In this appendix we present details of derivation of 
Eq. Pdj) which describes renormalization of g under non- 
equilibrium conditions in the weak-coupling regime. Ac- 
cording to general philosophy behind renormalization we 
successively integrate partition-function over high-energy 
components of field p. We split the scalar field into slow 
and fast components ip — > tp + where <p = ((p c , <p q ) and 
9 = (#c, 6q), and expand the action up to quadratic order 
in the fast field 9: 



S[p] -> S[p] 
6S 



bt[tp]6{t)dt+± I 9(t)K- t )[p}9{t')dtdt', 



h[<p] 



Stp{t) 



K: 



6 2 S 



6tp(t)8tp(t') 



)=0 



(Dl) 



Next we integrate out the fast components 9 and obtain 
the effective action for the slow components: 



S cS [p] = S[ip] - i J d ^yi b ii K "i"2 b "2 + |trlnJf 

(D2) 



= S[cp] - Si + S u . 



Here, frequencies u>i,ui2 lie in the energy window 
[A, A], A < A. The trace is understood to be over the 
frequencies in the same window as well as in the Keldysh 
space. High energy scale A in the AES-action is naturally 
set by the first term in Eq. (|T4| : A ~ gE c . Note that the 
linear in 9{t) term in (|D1[) does not generally disappear. 
But, as will be proven below, it is irrelevant since it leads 
to 1/A corrections. 

Next we perform the following decomposition 

K-^ip] = A' _1 [0] + [k~ 1 [<A - K- l [Q] 

= K-^O] +SI<- 1 [p] (D3) 

and treat the last term perturbatively. The operator 
A' _1 [0] determines a fast field propagator. It corresponds 
to perturbative Green function of the AES-action and fol- 
lows from Eqs. ([T6 ]) -([T? ]) : 

K R (t,t') = -i{ip c (t)ip q {t')}, K A {t,t') = -i{ip q (t)ip c (t')}, 
K K {t,t') = ~i{Vc{t)<Pc{t')). (D4) 

In the leading order the Wigner transform of the pertur- 
bative Green functions are given as 



K R = K At = - 
K K 



9 



2iImK R B u , 



(D5) 



where we neglect all time derivatives with respect to slow 
time since we are interested in high frequencies. The 
physical electron distribution function is bound to have 
sign-function as its limit at infinity F e — s- sgn(e), f-*oo. 
This yields the result 



K 



B, 



AQ 



2ni 



1 



K A ^ - 
v - w g uj + AQS 1 

E9 a Jde(l-F*F?_ u ) 



2g{uj + AQS) 



(D6) 



2 ^ 



de 



F" 



9o 



pa 



In general, AQ does not vanish. Next we find 
Sl = lJ^J dtdt'bitfK^e-^-^bit') 

dtdt'tfity-^-^K^Q] ( D? ) 



du)idu)2 



(2^)2 
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£ ' a exact Green's 
' function 



FIG. 7: The Dyson equation for pseudo-fermion self-energy. 



Performing integrations over fast frequencies u, u>i, 0J2 

we see that the first integral is ~ l/((i — t')A) and the 

2 

second one is ~ l/((£ — ti) (t' — t 2 )A ). Thus they are 
irrelevant for RG-analysis. It means that only term Su 
contains logarithmic in A corrections. 

As usual we are interested in the first non-vanishing 
(^-dependent correction: 



Sj!^ I tr V]}- (D8) 

Working out the trace in Eq. (|D8[) we obtain 



in 



dtdt' 



du 
2^ 



x (X c (t)X q (t) 



'A>|w|>A 



U R (t,t') 



U A (t,t') 
U K (t,t') 



X c (t>) 

x q (t') 



(D9) 



Substituting it into (|D2[) we see, that the structure of the 
AES-action is restored. The only difference is the change 
of the coupling constant given by Eq. (|41[) . Finally, we 
mention that in the case of non-zero AQ Eq. (|45)l should 
be changed to cj ~ maxje^, T r , TJ, |AQ|<5}. 



Appendix E: Renormalization of the pseudo-fermion 
action 



Here we provide details of the renormalization of the 
pseudo-fermion action (|rJ5|) which are used in Sec. I VII 



1. Renormalization of Z, A, and g. 

The exact pseudo-fermion Green's function can be 
written as 



1 



(El) 



Here, £ CT = Act/2 — ?y. To write it in the renormalized 
form (|75|) we redefine the theory's constants and write 
down the standard relations defining the Green's func- 
tion scaling Z , the renormalized gap A and the Green's 



function width Tjf respectively: 

Z=( 1 -<9 f Re£ R ' CT 



R,a\ 



igTl = -zZImEf'' 7 . 



(E2) 
(E3) 
(E4) 



To find the scaling Z and relate g, and A to their bare 
counterparts we solve the one-loop Dyson equation for 
the self-energy presented in Fig. [7j With the help of ((68)) 
we find 



£ 47T 

x Re 



dco 
2^ 



wB u Z(w) 



(E5) 



log — — iZ Im 



It is important to understand that scaling parameter Z 
cannot be put before the sign of an integral. Generally 
it is cut-off dependent and contains the factor ln(A/o;o), 
where A is an ultraviolet cut-off of the theory (E c in our 
case) while ujq is a characteristic scale of the Green's func- 
tion entering the integrand. To determine loq we notice 
that the integral in (|E5|) diverges, being determined by 
the behavior of the integrand in the large cj limit. That 
is why the characteristic scale of the Green's function 
entering (|E5|) is its running frequency: ujq ~ to. Solving 
Eqs. (|E5[) and (|E2|) with logarithmic accuracy we obtain 



1 dZ(A) _ g B(A)Z(A) 



Z 2 (A) dA At: 2 



(E6) 



A 



Integrating (|E6|) in the limits [u>o , E c ] we recover 
in complete analogy with equilibrium case. Thus, the 
renormalization procedure is outlined and the rest of for- 
mulae (l79l)-(l8Tl) are obtained in a similar fashion. 



2. Callan-Symanzik equation for (A/%/. 

The anomalous dimension 7 of (J\f) is introduced as 
^M P /(A,3) = Wp/(A,5,A). (E7) 

To extract 7 we write down the corresponding Callan- 
Symanzik equation for: (Af) p f(A,g,A) = J2 a (4'aip<y) ■ 
The tree-level (J\f) p f(A,g,A) is given by Eq. (ffl)). 
Following general strategy we write the corre- 
sponding Callan-Symanzik equation for the function 
(N) p f(A,g,A) in the form: 



P°Ig- + ^lk + ^' 



,<91nA 



(E8) 



where the corresponding ft - functions arc easily seen from 
Eq. (EU): 



4/ 



9 



P. 



A 



gA 

2tt 2 ' 



(E9) 
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FIG. 8: Correction to pseudo-fermion particle number {J\f) p f. 

The term with (3 g always contains extra g and can be 
dropped in the leading order. 

To find 7 we need to find (Af) p f in the next to the 
tree-level order. The diagram representing the correc- 



tion to pseudo-fermion particle number is presented in 
Fig. [5] Calculating with logarithmic accuracy and using 
extensively the fact that — s- sgn cj at large u> we obtain 

(AO p/ (A,g,A) = l- - F+ +- F ~ 

(E10) 

Plugging Eq. (|E"T0r) ) into Eq. (|E8]l) we find that 7 = 
that proves Eq. ((521) . 
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